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ABSTRACT 

Recent pulsar scintillation measurements from Ooty, in conjunction with those from Parkes and 
other large radio telescopes, are used for a systematic investigation of the Local Interstellar Medium 
(LISM) towards the general direction of the Loop I Bubble. For several pulsars, clear evidence is found 
for an enhanced level of scattering which is over and above what can be accounted for by the enhanced 
scattering model for the Local Bubble. These results are interpreted in terms of enhanced scattering due 
to turbulent plasma associated with the Loop I shell. Useful constraints are obtained for the scattering 
properties of the shell. The inferred value for the scattering measure for the Loop I shell is found to be 
~0.3 pc m~^°/'' . Assuming a shell thickness ^5-10 pc, this implies an average strength of scattering 
in the shell that is '^100-200 times larger than that in the ambient ISM. An alternative explanation, 
where the enhanced level of scattering is due to a possible "interaction zone" between the Local Bubble 
and the Loop I Bubble, is also considered; it is found to be somewhat less satisfactory in explaining 
the observations. The best fit value for the scattering measure for such an interaction zone region is 
estimated to be '^1.1 pc m^^^/-^ . 

Further, several pulsars beyond ~1 kpc are found to show enhanced levels of scattering over and 
above that expected from this "two-bubble model." For some of the low-latitude pulsars, this is found to 
be due to enhanced scattering from plasma inside the intervening Sagittarius spiral arm. We discuss the 
implications of our results for the interpretation of scintillation data and for the general understanding 
of the LISM. 

Subject headings: ISM:General - Structure - Bubbles - PulsarsiGeneral 



1. Introduction 

The Interstellar Medium (ISM) within a few hun- 
dred parsecs of the Solar System has been a topic of 
considerable observational and theoretical investiga- 
tion over recent years (e.g. Breitschwerdt, Freyberg 
& Trtimper 1998; Frisch 1996; Cox & Reynolds 1987). 
This region, often referred to as the Local Interstellar 
Medium (LISM), is known to contain several major 
features in the form of bubbles, supernova shells and 
Hi clouds. The Solar System itself is thought to re- 
side in a low-density. X-ray-emitting cavity of a mean 
radius ~ 100 pc, a region usually referred to as the 
Local Bubble (e.g. Cox & Reynolds 1987). It is be- 
lieved to be the remnant of a supernova explosion that 
occurred ~ 10*" yr ago. Prominent amongst other fea- 
tures in the LISM are the four large-diameter, almost 
circular rings seen in the all-sky distribution of radio 
continuum emission. These are appropriately named 



Radio Loops I to IV (Haslam, Khan & Meaburn 
1971) and are thought to be the projections of quasi- 
spherical bubbles. Of these, the Loop I Bubble — the 
largest in angular extent (solid angle ^ ^tt steradian) 
and also the brightest of the four loops — is of spe- 
cial interest (see Salter (1983) for a review). Its close 
proximity to the Sun allows an in-depth study; see 
for example, the work of Egger (1993) using data 
from the ROSAT soft X-ray Sky Survey, and that 
of Nishikida (1999) who combined the ROSAT PSPC 
data with IRAS Sky Survey and radio 21-cm data. In 
terms of origin. Loop I is thought to be an expand- 
ing super-bubble triggered by an epoch of star for- 
mation in the Sco-Cen association some ~ 10^ years 
ago. Further, it has been postulated that the proper- 
ties of the LISM may well be conditioned by the outer 
shock wave of this supernova remnant (Frisch 1981, 
1996; Lallement 1998). Besides these radio loops, the 
other well-known examples of nearby bubbles that 
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may be important in understanding the LISM are: 
the Eridanus bubble (distance from the Sun, D ^ 
100-150 pc), the Gum Nebula (D~200-250 pc), the 
Orion Bubble (D ~ 455 pc) and the Monogem ring 
(D ~ 100-1300 pc). 

Recent X-ray and UV data from ROSAT have led 
to several new insights into the structure of the LISM. 
A noteworthy result comes from ROSAT PSPC data 
suggesting an ongoing interaction between the Local 
Bubble and the Loop I Bubble. By comparing the 
diffuse X-ray background maps (in the 0.1 2.0 keV 
band) from the ROSAT All-Sky Survey (Triimper 
1983) and Hi data over a 160° x 160° region centered 
at {I, &)=(329°, +17.5°; the apparent centre of the Ra- 
dio Loop I), Egger & Aschenbach (1995; hereafter 
EA95; see also Egger 1993, 1998) recognize a "ring- 
like" structure between the Local Bubble and Loop I, 
where the Hi column density (Nh ) is even higher than 
that of the intervening dense Hi shell known to exist in 
the direction of the Sco-Cen OB association (Centu- 
rion & Vladilo 1991). From absorption- line studies of 
nearby stars, the distance to this interaction feature 
(where Njj jumps by nearly an order of magnitude) is 
estimated to be ~ 70 pc (see EA95), quite comparable 
to the distance to the neutral Hi wall (~ 40 ± 25 pc) 
inferred from ROSAT WFC star counts (Warwick et 
al. 1993) and optical and UV spectral line data (Cen- 
turion & Vladilo 1991). Interestingly, the formation 
of such a ring and the inferred density enhancement 
(by a factor ~ 25) are in good agreement with the 
numerical simulations of colliding interstellar bubbles 
(Yoshioka & Ikeuchi 1990). In this picture. Loop I 
is considered to be an active super-bubble, with at 
least one bubble having already formed a dense cool 
shell prior to collision. However, alternative inter- 
pretations do exist (e.g., Frisch 1996, 1998), wherein 
the Local Bubble is pictured as an appendix of Loop 
I. More data and modeling may help to distinguish 
clearly between the different scenarios. 

While much of our understanding of the struc- 
ture of the LISM comes primarily from X-ray and 
UV data, nearby pulsars are promising tools for en- 
hancing this understanding. Studies of dispersion and 
scattering of the radio signals from pulsars are useful 
means for probing the intervening ISM. Of these two, 
interstellar scintillation (ISS) effects are more likely 
to be influenced by the peculiar distribution of the 
ionized plasma along the line of sight (LOS) — such 
as clumps of enhanced density superposed on a uni- 
form distribution of ionized material. This is mainly 
because of (a) the nonlinear relation(s) of the scintil- 
lation properties to the electron density, and (b) the 
fact that scintillation effects depend critically on the 
relative location of the scatterer (or more generally. 



the actual distribution of scattering plasma along the 

sight line). For instance, if the medium is inhomo- 
geneous, in order to produce noticeable effects in the 
dispersion measure (DM), the density at the clumps 
has to be considerably larger than to produce similar 
effects in the scintillation data. However, if the pul- 
sar happens to lie within or near the clumped region, 
its effect on the scintillation properties could be sig- 
nificantly reduced. Thus, in some sense, dispersion 
and scintillation data can provide information com- 
plementary to each other. Nevertheless, (b) would 
imply that the interpretation of scintillation data is 
less immune to distance errors, and hence offers a 
better handle. 

It is quite plausible that the distribution of ionized 
plasma in and around large local features such as the 
Local Bubble and Loop I can considerably influence 
the dispersion and scintillation of nearby pulsars, and 
in some cases even scintillation of extra-galactic radio 
somccs. The large-scale distribution of free electrons 
and turbulent scattering plasma in the Galaxy has 
been modeled by Taylor & Cordes (1993; hereafter 
TC93). While sophisticated compared to its prede- 
cessors, the TC93 model takes very little account of 
the peculiar properties of the LISM due to individ- 
ual interstellar features. The recent years have seen 
an accumulation of observational evidence for the ef- 
fects of such features in the LISM. Early investiga- 
tions include the work of Phillips & Clegg (1992), who 
proposed that the scattering of radiation from the 
nearby pulsar PSR B0950-I-08 is probably dominated 
by weakly turbulent plasma present in the interior of 
the Local Bubble, and that of Hajivassiliou (1992), 
who invoked an ellipsoidal envelope of highly turbu- 
lent plasma to explain the directional anisotropy seen 
in the turbulent intensity maps derived from inter- 
planetary scintillation studies of radio sources. Even 
earlier, Rickard & Cronyn (1979) had suggested scat- 
tering from the outer shell of Loop I as the plausible 
cause of a statistically significant lack of interplane- 
tary scintillators seen in a band some 20° outside the 
brightest section of the Loop I radio emission, the 
North Polar Spur (NPS). 

In a previous paper (Bhat, Gupta & Rao 1998; 
hereafter BGR98), we presented a detailed, system- 
atic study of the LISM using pulsar scintillation data 
from the Ooty Radio Telescope (ORT). Our results 
and analysis strongly support the view that the scat- 
tering in the LISM is probably dominated by turbu- 
lent plasma at the boundaries of the Local Bubble. 
We proposed a simple model, wherein the Solar sys- 
tem is surrounded by an ellipsoidal shell-like struc- 
ture, with a size of ~100 pc in the Galactic plane and 
~500 pc in a plane perpendicular to this (Fig. 4). 
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The scattering structure has its center located at 
-20 35 pc from the Sun towards 215° < / < 240°, 
-20° < 6 < +20°. In this picture, the interior 
of the bubble is filled with plasma of relatively low 
turbulence (characterized by a scattering strength, 
~ lO"** m~^°/^), whereas the shell material 
(thickness ~1-10 pc) has a scattering strength —50- 
800 times larger than that in the ambient ISM (in- 
tegrated strength of scattering, scattering measure, 
given by O.IK SMlb <0.28 pc m-^o/a ). The contri- 
bution of the shell thus dominates the total scatter- 
ing, which would imply that the scattering geometry 
towards many sight lines can be approximated by a 
"thin screen" placed at the bubble boimdary (location 
in the range —20-200 pc). This model successfully 
explained the enhanced level of scattering measured 
towards a number of nearby pulsars. 

The Ooty experiment covered only a few pulsars 
that would be useful for a study of Loop I. How- 
ever, the recent results from Parkes observations of 
a large number of southern pulsars (Johnston, Nicas- 
tro & Koribalski 1998) has significantly improved the 
ISS data available for probing the ISM in and around 
Loop I. This motivated us to take a more detailed 
look at the distribution of scattering material towards 
Loop I and beyond, a study which forms the main 
theme of this paper. The paper is organized as fol- 
lows: § 2 describes observational data and our anal- 
ysis techniques; § 2.4 and § 2.5 present the modeling 
of the scattering plasma associated with Loop I. The 
uncertainties relevant to our analysis are discussed in 
§ 3.1. In § 3.2, we consider some possible alternative 
models, while later sections of § 3 discuss some gen- 
eral implications of our results for scintillation data 
and for the LISM. In § 4, we summarize our conclu- 
sions. 

2. Observational Data, Analysis Techniques 
and Modeling 

2.1. Sample Selection 

For the present analysis, we are interested in pul- 
sars whose scintillation properties are likely to be in- 
fluenced by the structure of Loop I. In order to pre- 
select pulsars useful for this purpose, we adopt a ge- 
ometry and size for Loop I from the published lit- 
erature. Specifically, Berkhuijsen, Haslam & Salter 
(1971) derived an angular diameter of 116° ± 4°, and 
a centre at {I, b) = (329.0° ± 1.5°, -hl7.5° ± 3.0°; see 
also EA95). Further, on the basis of the plausible con- 
nection between the origin of Loop I and the Sco-Cen 
OB association (e.g. Egger & Aschenbach 1995), it is 
fair to assume the loop center to be near the center 
of mass of the association (—170 pc). Based on the 



above, we model the Loop I Bubble as a spherical shell 

of size —290 pc, with the center located at —170 pc 
towards (329°, 17.5°). This yields 52 pulsars within 
— 2 kpc of the Sun whose sight lines intersect the pro- 
jected area of Loop I on the sky. On carrying out a 
literature search, we find scintillation measurements 
to be available for only 20 of these (see Table 1). The 
majority of these scintillation measurements are from 
recent observations with the Parkes and Ooty radio 
telescopes, while a few are from Arecibo (data from 
Johnston et al. 1998; Bhat et al. 1999b, Gothoskar & 
Gupta 2000, and Cordes 1986). The measurements of 
decorrelation bandwidth [vd^meas ) and the observing 
frequencies [ fobs ) are listed in columns (5) and (6) of 
Table 1, respectively. The distance estimates in col- 
umn (4) are derived from dispersion measures (DMs) 
and the model of Galactic electron density (rie ) by 
TC93, except for PSRs J1456~6843 and J1744-1134. 
For these two pulsars, we use independent distance 
estimates derived from the measurement of annual 
trigonometric parallax (Bailes et al. 1990; Toscano et 
al. 1999). 

Although the decorrelation bandwidth measure- 
ments listed in Table 1 are at different observing fre- 
quencies, we have scaled them to a common frequency 
of 327 MHz, assuming a Kolmogorov scaling law (i.e., 
a wavenumber spectrum with a slope of 11/3 over 
the spatial scales of interest; Vd oc. frequency*'^). We 
note that this can potentially produce some errors, as 
the exact nature of the electron density wavenumber 
spectrum is still debated. However, there is substan- 
tial observational evidence in favor of an a « 11/3 
spectrum towards many sight lines in the LISM (e.g. 
Bhat, Gupta & Rao 1999a). Even if this were not 
strictly correct, most measurements from Parkes will 
be only marginally affected by an incorrect frequency 
scaling. This scaling bias may be significant for 
Vd values at 1.5 GHz; however, these are for relatively 
distant objects (D>1.2 kpc), and hence less critical 
for the investigation of scattering in the LISM. 

2.2. Distribution of Scattering: Choice of the 
Method 

There are three different methods by which pul- 
sar scintillation measurements can be used to in- 
vestigate the distribution of scattering material in 
the Galaxy: (1) using the decorrelation bandwidth, 
Vd (or its equivalent, the temporal pulse broadening 
time, Tp ) — this quantifies the scattering measure 
(SM) which characterizes the total amount of scatter- 
ing along the line of sight. This technique has been 
used extensively to investigate the large-scale distri- 
bution of Cn in the Galaxy (Cordes, Weisberg & Bo- 



3 



riakoff 1985; CWB85 hereafter; Cordes et al. 1991; 

TC93): (2) using measurements of angular broaden- 
ing {Oscatt ) in conjunction with Tp — wherein the 
differing weighting functions of the two observables 
can be used to determine a more exact distribution 
of Cn along the LOS (e.g. Gwinn, Bartel & Cordes 
1993): (3) using the hybrid method recently proposed 
by Cordes & Riekett (1998) wherein the diffractive 
scintillation measurements (decorrelation bandwidth 
and scintillation timescale, ), in conjunction with 
the pulsar proper motion and distance, can be used 
to obtain the distribution of C^ along the LOS (e.g. 
Chatterjee et al. 2000). 

Of these methods, (2) is not relevant in our case, as 
measurements of Oscatt and Tp do not exist for bulk 
of the objects in Table 1. To date, proper motion 
measurements arc known for 9 of the 20 pulsars in 
Table 1. However, uncertainties are too large for the 
values to be useful in most cases, and we were unable 
to derive any meaningful results on the distribution 
of scattering in and around Loop I. We therefore re- 
stricted ourselves to method (1) for the analysis de- 
scribed in this paper (see Appendix A for a detailed 
description of the method). 

2.3. Comparison with the Local Bubble 

model 

We first examine how well the scintillation data in 
Table 1 agree with the predictions of the Local Bubble 

model of BGR98, as described in § 1. Fig. 1 shows a 
plot of the ratio of the measured to predicted decorre- 
lation bandwidths {i'd,meas/^d,pred) against the pul- 
sar distance estimates. We introduce a qiiantity, eaii , 
that is a measure of the degree of agreement between 
'^d,meas and i'd,pred • This is expressed as 

1 r / \ 1 2 

eall = ^ r 2^ log -f (1) 

where Np is the total number of pulsars for which the 
comparison is being made. The logarithm has been 
taken to give equal weight to discrepancies that are 
below and above unity when computing eaii ■ ^ 

As is obvious from Fig. 1, there are large discrep- 
ancies (ranging from a factor 2 to as much as ^50) 
for the case of the Local Bubble model (shown by 
crosses). Considering that most measurements in Ta- 
ble 1 are not from observations averaged over many 

We will use the quantities e<ifcpc a-nd f>ifcpc when referring to 
those parts of the data that correspond to objects with D S 1 
kpc and D <; 1 kpc, respectively. 



epochs, a potential explanation for part of these dis- 
crepancies may be errors due to long-term refractive 
interstellar scintillation (RISS) effects (e.g. Gupta, 
Riekett & Lyne 1994; Bhat et al. 1999a). However, 
most discrepancies arc considerably larger than the 
worst case factor of 3-5 that can be accounted for by 
RISS effects. Moreover, it is striking that all the ra- 
tios are less than unity, implying that the scattering 
strengths are systematically larger than that can be 
accounted for by the Local Bubble model. Further, 
a systematic trend with distance is also evident in 
Fig. 1, whereby i'd,meas / '^d,pred (for the Local Bub- 
ble model) tends to lie in the range 0.01-0.1 for pul- 
sars beyond 1.2 kpc, but between 0.1 and 1 for those 
nearer than this. In contrast, for a sample of pulsars 
within 1 kpc of the Sun in the complementary sky 
(shown by □ in Fig. 1, cf. BGR98), most ratios fall 
within a factor 2-3 of unity. The value of eaii for this 
sample (« 0.07) is about 16 times lower than that for 
the Local Bubble model applied to the current data 
set (caii ~ 1.1, as shown in Table 3). The simplest 
interpretation is that there is yet another source of ex- 
cess scattering (in addition to the Local Bubble shell) 
in the general direction of Loop I. 

2.4. Enhanced Scattering from the Shell of 
the Loop I Bubble? 

A number of previous studies have revealed the 
existence of strong excess scattering towards several 
sight lines in the Galaxy (e.g. Moran et al. 1990). In 
the context of modeling the large-scale distribution 
of C^ in the Galaxy, Cordes et al. (1991) propose 
a 'clump' component, possibly associated with super- 
nova shocks or Hii regions, to explain such unusually 
high scattering. These are regions of intense turbu- 
lence, with the strength of scattering many orders of 
magnitude larger than that in the typical ISM. How- 
ever, the mean free paths expected for such clumps 
are rather large (typically ^5 10 kpc). There seems 
to be no observational evidence for the presence of any 
such regions in the LISM. Strong excess scattering 
may also arise if the sight line to a pulsar intercepts 
an Hii region (e.g. Gupta et al. 1994). To the best 
of our knowledge, none of the pulsars in Table 1 rep- 
resents such a situation. Even if the LISM contains 
strong scattering regions that are hitherto unknown, 
these could be relevant only for a few lines-of-sight, 
whereas the data in Fig. 1 suggest a common source 
of enhanced scattering for many sight lines. Given the 
current understanding of the LISM, and also on the 
basis of the results from our earlier work, we postu- 
late that this enhanced scattering is probably caused 
by the shell of the Loop I Bubble. 
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2.4- 1- Modeling the Scattering due to the Loop I 
Shell 

Here we try to explain the observed discrepancies 
in decorrelation bandwidths for nearby pulsars in the 
direction of Loop I in terms of enhanced scattering 
from a combination of the Local Bubble and the Loop 
I Bubble. A full modeling of the problem involves the 
determination of the parameters characterizing their 
size and location, and the distribution of scattering 
in and around them. This was done by BGR98 for 
the case of the Local Bubble alone. However, in the 
present case, since we already have a model for the 
scattering from the Local Bubble and there are very 
good models for the physical dimensions of Loop I, 
we do not really need to do such full-blown modeling. 
Instead, we take the following simpler approach. 

First, we assume a scattering model for Loop I 
which is very similar to the multi-component model 
for the Local Bubble, i.e., an interior filled with 
a plasma of relatively low turbulence (the LOS- 
averaged strength of scattering. C^^ ^ 10^^'^ -j 
and a dense, highly turbulent shell in which the 
scattering strength is many times larger than in the 
ambient ISM that lies beyond Loop I. For this am- 
bient medium, we assume ~ 10~'^-^ m~^°/^ (see 
BGR98). 

Second, as described in § 2.1, the size and structure 
of Loop I are fairly well constrained by previous work 
(e.g. Berkhuijsen et al. 1971). We simply adopt these 
values for characterizing the physical dimensions and 
location of the Loop I shell in our analysis. 

Third, we simplify the problem by assuming the 
thickness of the Loop I shell {dLi ) to be small 
compared with the pulsar distances involved. The 
problem then essentially reduces to obtaining the 
best estimate of the scattering measure due to the 

Loop I sheU ( SMli = J^" Cl^Lish (0 dl, where 
i/sh denotes th(> scattering strength inside the 
shell) that will minimize the discrepancies between 
i>d,meas and i^d,pred ■ Notc that the spherical shell ge- 
ometry of Loop I implies that SMli should be a 
function of the LOS, or more precisely, of the angular 
distance of the LOS from the loop centre (9). This 
is because the actual path length through the shell 
intercepted by the LOS varies with O. The effect will 
be more prominent for sight lines that are close to 
being "tangential." In the analysis that follows, we 
will constrain SM as defined for a LOS that is normal 
to the shell (i.e., along = 0). 

Fourth, instead of obtaining a solution that will 
give a best case minimum for all the pulsars in Ta- 
ble 1, we prefer to solve for the SM of the Loop I shell 



using only two pulsars PSRs J1744-1134 and J1456- 

6843. This is justified on the grounds that (a) for 
many pulsars in Table 1, the distance uncertainties 
and the large errors in the estimates of the decorre- 
lation bandwidth (due to RISS effects), make it less 
meaningful to do a global fit, and (b) for several pul- 
sars at distances z 1 kpc. Loop I may not necessarily 
be the only source of enhanced scattering and a global 
fit for SMli could therefore be perturbed by this ef- 
fect. 

The choice of the above two pulsars is dictated by 
the following arguments. Both pulsars are heavily 
scattered (Johnston et al. 1998) and have precise, in- 
dependent distance estimates from parallax measure- 
ments — 357l3^ pc for PSR J1744-1134 (Toscano et 
al. 1999) and 4551^^ pc for PSR J1456-6843 (Bailes 
et al. 1990). Interestingly, the estimate for PSR 
J1744-1134 is over twice the value of 166 pc derived 
from the model of TC93. The independent distance 
estimates allow the scattering geometry along these 
lines of sight to be fairly well constrained. This is 
illustrated in Fig. 2, where the expected locations 
of the Local Bubble (Db,LB ) and Loop I boundaries 
(Dbi.Li and Db2,Li) in these directions are clearly 
marked. The farther boundary of the Loop I shell 
(L)b2Xi) lies very close to the mid-point to both ob- 
jects, and is thus optimally placed for maximum con- 
tribution to the observed decorrelation bandwidth. 
Such a geometry will be consistent with the enhanced 
level of scattering inferred for these two pulsars. 

Appendix A describes our method for comput- 
ing the expected decorrelation bandwidth {i'd,pred ) 
for a given distribution of the scattering material 
along the LOS. We note that if homogeneously dis- 
tributed scattering material (i.e., equation [Al], with 
the integral replaced by C^D, where denotes the 
LOS-averaged strength of scattering) were to ex- 
plain the observed scattering for the lines of sight of 
these pulsars, the implied level of scattering strength 
(C2=10-2-« m-20/3 towards PSR J1744 1134, ob- 
tained from the value reported by Johnston et al. 
[1998], after scaling for the new distance) would be 
several times larger than that expected for the typical 
ambient ISM. This seems quite unlikely. The LISM 
as relevant to the Local Bubble model of BGR98 can 
be treated as an inhomogeneous scattering medium, 
which can be best represented by multiple compo- 
nents of different strengths of scattering ( ) along 
the LOS (see equation [A2], for example). This can be 
represented by a simpler, modified version of c;quation 
(A2) (i.e., the right hand side consisting only the first 
two integrals and the last integral with the lower limit, 
Db2,Li+ dLi , replaced with Db,LB + c^lb )• If we were 
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to attribute the entire excess scattering to material in- 
side the Local Bubble shell (located at ^30-60 pc in 
this direction; see Fig. 2), then C^j^^^^ would have 
to be over an order of magnitude larger than that de- 
rived by BGR98; alternatively, the shell would have 
to be much more extended in this direction (say, by 
an order of magnitude). Again, either of these possi- 
bilities is rather unlikely, as it would require a drastic 
change in the properties of the Local Bubble in this 
direction. Hence the possibility that the Loop I shell 
is the source of the enhanced scattering is most plau- 
sible. 

We now estimate the optimal value of SMli by 
modeling the scattering towards PSRs J1744 1134 
and J1456-6843, taking into consideration the scat- 
tering due to the Local and Loop I bubbles, as well 
as that due to the distributed plasma (i.e., ambient 
ISM) beyond Loop L Our model also takes into con- 
sideration the characteristic fall-off with z-height for 
the ambient scattering plasma (scale height ^^500 pc; 
see TC93), and a Gaussian fall-off (scale height ~125 
pc; see BGR98 for arguments) for the shell material. 
The interiors of both bubbles are taken to be filled 
with weakly turbulent plasma (C^ ~ lO"''-^ m^^°/^ ). 
The shell of the Local Bubble is taken to have a 
scattering measure ( SMlb ) of 0.2 pc m"^"/"^ (i.e., 

~ a few 100 times larger than that in the ambient 
medium). The method for computing the expected 
decorrelation bandwidth [vd.pred ) for such a geom- 
etry is described by equation (A2) in Appendix A. 
To determine SMli , we start with an initial value 
of 0.2 pc m^^"/^ , and vary this from 0.02 to 2.0 
pc mT^^/^ in steps of 0.01 pc m~^°/^ . Interestingly, 
the best agreement between the predicted Vd values 
and those measured is obtained for SMli ~ 0.29 
pc m"^"/"^ , a value quite comparable to that inferred 
for the Local Bubble shell. 

Using the SM of the Loop I shell, as constrained 
by the above method, we compute new values for the 
Vd ratios. These are shown by the filled star sym- 
bols (★) in Fig. 1. It is quite remarkable that the 
above modeling successfully removes the discrepan- 
cies for pulsars out to a distance of a few 100 pc. The 
Call value for this model (« 0.55) is a significant im- 
provement over that for the earlier model of the Local 
Bubble alone, as listed in Table 3. Further, as can be 
seen from this table, there is a dramatic improvement 
in e<ifepc (i-e., for pulsars within 1 kpc), and a much 
smaller one in e^ikpc- exactly as expected for the 
Local Bubble -|- Loop I Bubble model. 



2.5. Enhanced Scattering from the Sagittar- 
ius Spiral Arm 

Next we address the systematic downward trend 
for the ratio fd,meas/t'd,prerfin Fig- 1 for pulsars be- 
yond ~1.2 kpc. A closer examination of these LOSs 
suggests material within the Sagittarius spiral arm to 
be the most plausible source of enhanced scattering 
for several of the pulsars. We show here that these 
results can be better explained by taking into account 
the enhanced level of scattering that can be expected 
due to this spiral arm as per the existing model of 
TC93. With the arm parameters of the TC93 model 
and the distance estimates based on this model, the 
sight lines of six of our pulsars pass through this spi- 
ral arm. To model the effect of this, we have simply 
adopted the arm locations and density parameters of 
TC93. They consider a squared hyperbolic secant for 
the z-depcndence and a Gaussian fall-off for the radial 
dependence for the arm density, with 300 ± 100 pc for 
the scale height and 300 pc for the half width. With 
an arm density of 0.08 cm^"^ i.e., several times larger 
than the typical (rie) in the LISM, and a fluctuation 
parameter (F) of 6^2) the implied level of scatter- 
ing is some 2 orders of magnitude larger inside the 
arm ( C^^^^ ~0.05 m'^o/^ ) than in the LISM. Hence, 
substantial amounts of enhanced scattering are ex- 
pected for these pulsars. On incorporating this, the 
discrepancies are further reduced (see the open circles 
in Fig. 1), with a significant improvement in overall 
agreement (eaii « 0.22) as well as for e^ikpc (see Ta- 
ble 3). The best agreement is seen with the lower 
and upper limits allowed (by the model of TC93) for 
the values of the fluctuation parameter and the scale 
height, respectively. 

3. Results cmd Discussions 

We have shown that the combined effect of en- 
hanced scattering from the Local Bubble shell, the 
Loop I Bubble shell and the Sagittarius arm goes a 
long way towards explaining the observed scattering 
properties of pulsars in the general direction of Loop 
I. Our new model consists of: (a) an ellipsoidal shell of 
SM~0. 1-0.3 pc m^^"/'^ to account for scattering due 
to the Local Bubble, (b) a spherical shell of SM~0.3 
pc m~^°/^ to characterize the scattering due to Loop 
I, (c) the ambient ISM (Cf- lO^^-S m-20/3 ) the 
inter-arm region, and (d) the Sagittarius spiral arm 
(C^-^ 0.047 m~^°/^). We now discuss some of the 
implications of these results. 
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3.1. Scattering from the Loop I Shell 

The agreement between the measurements and the 
predictions of our new model is such that the scatter- 
ing properties of most pulsars out to a distance of 1 
kpc are successfully explained (fd,meos/i'd,predin the 
range 0.5-2.0 for 8 out of 10 objects). The major 
outliers are PSRs J1713+0747 and J1730-2304, with 
residual discrepancies of ~ 4 and ~ 3 times respec- 
tively between the predicted and measured values (see 
Fig. 1) that remain to be explained. These discrep- 
ancies are possibly due to RISS, or distance errors, 
or perhaps as yet im- modeled source (s) of enhanced 
scattering. For more distant pulsars (D S 1.2 kpc), 
although the discrepancies are better accounted for, 
there still seems to be a systematic downward trend 
with distance, which is accounted for by extra en- 
hanced scattering due to the Sagittarius arm. The 
results thus clearly substantiate the role of Loop I as 
a source of enhanced scattering in the LISM. 

We now address the various sources of uncertain- 
ties relevant to our analysis. These include: (a) the 
uncertainties in size and geometry of (i) the Local 
Bubble, and (ii) Loop I, (b) the uncertainty in the 
(integrated) scattering strength of the Local Bubble 
( SMlb), (c) distance uncertainties, and (d) errors 
on the measurements of Ud . First, we note that the 
Local Bubble geometry itself is not very well deter- 
mined (see BGR98), and the expected location of the 
boundary is in the range ~24-58 pc for PSR J1744- 
1134, and -44-74 pc towards PSR J1456 6843 (see 
Fig. 4). Taking this into consideration, SMli could 
be in the range 0.25-0.3 pc m~^°/^ . However, the 
uncertainties in the angular size and location of Loop 
I itself are small (see Berkhuijsen et al. 1971), hence 
we do not expect them to affect SMli apprecia- 
bly. Turning to (b), the strength of scattering of 
the Local Bubble shell ( SMlb ) itself has a signifi- 
cant uncertainty : O.IK SMlb <0.28 pc m-^°/^ (see 
BGR98). However, this seems to affect SMli only 
marginally (0.24-0.32 pc m-^o/S). The distances to 
PSRs J1744-1134 and J1456-6843 have ~10% uncer- 
tainties, which translate in to a variation of 0.28-0.33 
pc m~^°/'^ for SMli . As for the measurements of 
Vd , the dominant source of errors are due to their 
apparent variations caused by RISS effects on time 
scales from days to weeks. However, these uncertain- 
ties are hard to estimate. Nevertheless, if we consider 
a typical case of a factor of two variation in the mea- 
sured values of i/^ , SMli may lie in the range 0.15- 
0.57 pc m~^°/^ . Hence, the RISS-induced errors in 
measured deeorrelation bandwidths are likely to be 
the dominant source of errors in the final estimate of 
the strength of scattering of the Loop I shell material. 



3.2. Alternative scenario: Scattering from an 
Interaction Zone between the Local Bub- 
ble and Loop I? 

In the above analysis, we have attempted to un- 
derstand the observations in terms of enhanced scat- 
tering caused by turbulent plasma within the shells 
of the two bubbles. Although this two-bubble model 
seems to successfully explain the scintillation mea- 
surements for most nearby pulsars, there is an in- 
teresting alternative that deserves consideration. As 
described in § 1, the work of EA95 and Egger (1998) 
suggest that the two bubbles are probably undergo- 
ing a collision process; the observational evidence for 
a dense interaction feature at ~70 pc towards Sco- 
Cen (Fig. 4 of EA95) supports this. Interestingly, 
the expected location of the closer boundary of Loop 
I (Dbi,Li) is either within the range of or near the 
Local Bubble boundary (Db,LB ) as constrained by 
BGR98 (see Fig. 4). In particular, for several lines- 
of-sight, Dbi,Li lies in between the smaller and larger 
boundaries (Dj^ and Dj^ respectively) of the 
Local Bubble. The possibility of enhanced scattering 
from an "interaction zone" between the two bubbles 
is therefore worthy of consideration. 

To examine this possibility, we consider a model 
where the measured level of enhanced scattering 
is entirely attributable to an "interaction wall" or 
"zone" whose effective location is at the mean dis- 
tance of the LB and Loop I boundaries. The 
two extreme possible geometries for the LB would 
thus imply two possible locations for the inter- 
action zone (IZ): =0.5(D^ ^b +Dbi,Li) and 

j2 =0.5(D^ LB +l^bi,Li), which correspond to the 
cases of nearer and farther boundaries of the LB, re- 
spectively. These are listed in column (5) of Table 2 
(db,iz =Db,iz /D). In the discussion below, we will 
refer to these as the IZ-A and IZ-B models, respec- 
tively. We also assume that the thickness of this zone 
{diz ) is comparable to that of the LB or LI shell, and 
also that it is much smaller than the pulsar distances. 
Further, like the case of the Loop I shell, the path 
length of the sight line through a possible interac- 
tion zone region, and consequently its contribution 
to the scattering measure, will be a function of the 
angular distance of the LOS from the loop centre. In 
order to determine the scattering measure of the IZ 
( SMiz ) , we performed a rather similar analysis to 
that described in § 2.4.1, in which PSRs J1744-1134 
and J1456-6843 are used to estimate the value for 
SMiz . For the IZ-A model (Dbjz =Db.iz the best 
agreement between the measured and predicted val- 
ues for Vd is obtained for SMiz «1.12 pc m~^°/^ , 
which is -^4 times larger than the value of SM for 
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the Loop I shell, and ~4-10 times larger than that 

for the LB. The required SM for the IZ-B case is 
somewhat smaller: ~3 times the value for SMli , 
and ~3-9 times SMlb • This is mainly due to the 
relatively further location of the zone in this case. 
Using these values for SMjz , we re-computed new 
values for Vd,meas/t^d,pred ■ The results for the IZ-A 
model are shown in Fig. 3. The agreement is not as 
good as that achieved with the two-bubble model, as 
indicated by somewhat larger values for the three e 
parameters that quantify the level of agreement (see 
Table 3). On a closer comparison of Figs. 1 and 3, 
it is obvious that the agreement is somewhat poorer 
for several objects at D 5 1 kpc. Nevertheless, the 
IZ model(s) may still be considered as a possible al- 
ternative to the two-bubble model. However, unless 
there are well studied objects located in the Loop 
I interior, it is not easy to distinguish conclusively 
between these two scenarios. 

As a natural extension of the IZ model(s), we in- 
corporate scattering due to the Sagittarius spiral arm 
(as discussed in § 2.5) in order to improve upon the 
agreement for pulsars beyond 1 kpc with sight lines 
intercepting that arm. The results for the IZ-A-t-SA 
case are shown in Fig. 3. The overall agreement, as 
indicated by eaii =0.27 and 0.25 respectively for the 
IZ-A-t-SA and IZ-B-I-SA cases, is somewhat poorer 
than that achieved with the two-bubble+SA scenario 
{eaii =0.22, see Table 3). 

3.3. Discussion 

3.3.1. Implications of the Results 

(a) Interpretation of the Scintillation Data 

From the results of this work it is quite evident that 
the structure of the LISM plays an important role in 
the interpretation of the scintillation data. The outer 
shells of the Local Bubble and Loop I may very well 
be the dominant sources of scattering towards many 
directions in the LISM, especially at higher Galactic 
latitudes. Interestingly, independent evidence in favor 
of similar close-by scattering screens (~25-250 pc) 
comes from recent observations of cm-wave ISS and 
intra-day variability (IDV) of quasars (e.g., Dennett- 
Thorpe & de Bruyn 2000; Rickett 2000). In partic- 
ular, the screen location of ^^25 pc and strength of 
scattering, RiO.2 m"^"/'^ , inferred by Dennett- 
Thorpe & de Bruyn (2000) towards the IDV quasar 
J1819+3845 is in very good agreement with the lo- 
cation of the Local Bubble boundary (~25-50 pc) 
and the scattering strength ( SMlb ) expected in this 
direction. Further, from the work of Hjellming & 
Narayan (1986), the bulk of the scintillation (RISS) 



of the radio source PKS 1741-038 at 1.49 GHz is 

caused by a single screen located at ^^140 pc from 
the Sun. Interestingly, the LOS to this object (/ = 
21°. 6, 6 = 13°. 1) lies very close to the "inner ridge 
feature" of the NFS (beginning at I « 22°, 6 = 14°). 
Also, Hjellming & Narayan infer that there is signif- 
icant excess scattering ( ~ 10~^'^ m~'^^/^ ) along 
this LOS. Recently, Lazio et al. (2000) report multi- 
epoch VLBI observations of this object as it under- 
went an extreme scattering event (ESE); the inferred 
level of scattering for the ESE lens is orders of magni- 
tude larger than that in the ambient medium. Given 
the close proximity of the object to Loop I, it is quite 
probable that the structure (s) that caused the ESE 
are located at the Loop I shell. Similarly, if an inter- 
action wall between the two bubbles exists, it could 
potentially act as a thin scattering screen for nearby 
pulsars with sight lines within 270° < I < 30° and 
—40° < 6 < -1-80°. Such close-by scattering screens 
may also explain the shorter than expected timescales 
seen with the slow intensity variations of pulsars and 
the low frequency variables (Gupta, Rickett & Coles 
1993; Spangler et al. 1993). In addition, the enhanced 
scattering from such bubble shells may be responsi- 
ble for some of the unusual scattering effects. In this 
context, it is interesting to note that the location of 
the scatterer that caused the multiple imaging event 
of PSR B1133-M6 in the Ooty data (Gupta, Bhat & 
Rao 1999) was found to match well with the expected 
location of the Local Bubble shell (Db,LB ~ 0.77 D) 
in this direction. Similarly, it has also been suggested 
that the the structures that cause extreme scattering 
events in the radio light curves of some quasars are 
probably associated with the Galactic loops (Fiedler 
et al. 1994). Recent work by Toscano et al. (1999) 
and Chatterjee et al. (2000) provides some evidence 
for enhanced scattering in the third galactic quadrant 
probably associated with an interface region between 
the Local Bubble and the GSH 238-^00-^09 super- 
bubble. All of the above arguments clearly support 
the view that the structure of the LISM needs to be 
taken into account for a proper interpretation of dis- 
persion and scintillation data. 

(b) On the Nature of the Bubble Shells 

In our earlier work (BGR98), we argued that the 
elongated ellipsoidal shell of turbulence derived from 
the Ooty data is possibly associated with the Local 
Bubble. The analysis described in § 2.4 presents quite 
convincing cnddence for an outer turbulence shell for 
Loop I also (the existence of a turbulent interaction 
zone is an alternative possibility, though). Thus, from 
observational data, it appears that interstellar bub- 
bles in general may have turbulent outer shells. 
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On the theoretical front, a number of authors 

have dealt with the evolution of bubbles in the ISM 
(Ikeuchi 1998 and references therein); in particular, 
special attention has been paid to the major local 
features such as the Local and Loop I Bubbles. Inter- 
estingly, it turns out that turbulent outer shells can 
indeed be expected for interstellar bubbles. Specifi- 
cally, the recent model of Breitschwcrdt, Freybcrg & 
Egger (2000) predicts turbulence to exist at the in- 
terface region of the Local and Loop I Bubbles as a 
natural consequence of the hydro-magnetic Rayleigh- 
Taylor instability caused by the interaction of the two 
bubbles. A similar mechanism may also apply for the 
Loop I shell. In an earlier work, Breitschwcrdt & 
Kahn (1989; see also Kahn & Breitschwcrdt 1988) 
showed that turbulent shells could also be expected 
for stellar-wind-blown bubbles as a consequence of 
acoustic instability. So, in general, the shells around 
bubbles may be expected to be turbulent. There are 
no explicit predictions available for the level of turbu- 
lent intensity expected at the shell or in interaction 
zone regions, or for a relation connecting the turbu- 
lent intensity to the physical parameters governing 
the process(es). Nevertheless, the basic results from 
our investigations using pulsar data appear to be in 
accordance with some of the theoretical models of 
bubble evolution. 

3.3.2. Estimation of the Shell Density 

The strength of scattering we derived for the outer 
shell of Loop I is substantially higher (^100- 500 times 
depending on the shell thickness) than that measured 
towards typical LOSs in the LISM. The regions of en- 
hanced electron density (rie ) arc also thought to be 
plausible sites for enhanced scattering (e.g. the Gum 
Nebula) , though the exact relationship between these 
two is not clearly known. If we consider to be 
uniform within the shell or the interaction zone re- 
gions (i.e., SMli .= C^xish d-Li , SMjz = , 
where (Ili and diz represent the path lengths through 
the shell and interaction zone regions respectively), 
and assume a simple relation ( oc nl ), then an in- 
direct estimate of electron density in the shell (ne,s/i ) 
is possible for a given shell thickness. Considering the 
Loop I shell (SM-G.3 pc m-^o/s ), for a shell thick- 
ness ~1-10 pc, this yields a density ~10-30 times 
larger than the ambient ISM value. Somewhat higher 
values (~20-50(ne) ) result for the interaction zone 
models described in § 3.2. Interestingly, a similar level 
of density enhancement is seen for the neutral gas at 
the annular interface region of the two bubbles (EA95; 
Egger 1998). For the above densities, the effects on 
dispersion due to the shell region will be significant 
only for pulsars within a distance of a few 100 pc. 



Taking this into consideration would, however, result 

in a much lower ambient density ((ng) «0.GG7 cm~^ ) 
towards PSR J1744-1134, the disk pulsar with the 
lowest known Ue . For a similar density enhancement 
(~10 times) for the shell region towards PSR J1456 - 
6843 (the only pulsar in quadrant 4 with a measured 
parallax), consideration of shell dispersion would re- 
sult in an ambient density of « 0.016 cm~^ for the 
two-bubble model, and ~ 0.013 cm""^ for the in- 
teraction zone model. These values are comparable 
to (ue) measured towards several objects in quadrant 
1 with known parallaxes (see Toscano et al. 1999). 
Therefore, it is quite plausible that the Loop I shell 
(or the interaction zone) contributes significantly to 
the dispersion of nearby pulsars in quadrant 4. 

3.3.3. Electron Densities towards PSRs J 1744^1 134 
and J 1456-6843: Evidence for a Dense Wall? 

It is interesting to note that for the two pulsars 
- PSRs J1456-6843 and J1744-1134 - with indepen- 
dent distance estimates, the mean electron densities 
along the sight lines differ significantly: (rie) towards 
the former is almost twice that for the latter, despite 
their comparable distances and z-heights (67 and 57 
pc respectively). A closer look at their LOSs vis-a- 
vis the geometries of the Local and Loop I Bubbles 
(Figs. 2 and 4) indicates the following: (1) the Local 
Bubble cavity extends out almost twice as far towards 
PSR J1456-6843 as towards PSR J1744 1134 (~44- 
74 pc and ^24 58 pc respectively, for the geometry 
in Fig. 4), (2) the interior of Loop I covers a much 
longer section of the LOS towards PSR J1456-6843 
(Dint=Db2.Li-Dbi.Li ~ 233 pc for J1456 6843, com- 
pared to 157 pc for PSR J1744 1134; see Fig. 2), and 
(3) the extent of the ambient medium beyond Loop I 
is not much larger along the sight line towards PSR 
J1456 6843 (-192 pc compared to -159 pc for PSR 
J1744-1134). All of these make it harder to explain 
the larger (ng) towards PSR J1456 6843. 

There appears to be two plausible explanations for 
the above: (a) the ambient ISM in quadrant 4 is sim- 
ply much denser than that in quadrant 1. However, 
in the absence of other pulsars in quadrant 4 with 
measured parallax, there seems to be no independent 
way to confirm this, (b) The existence of a dense 
interface region between the Local and Loop I Bub- 
bles, as described in § 3.2. On a closer inspection of 
the sight lines of the two pulsars, case (b) seems the 
more likely. The LOSs towards PSR J1744-1134 in- 
tercepts the "annular Hi ring," as recognized by EA95 
(sec also Egger 1998), whereas the LOS toward PSR 
J 1456-6843 lies well within the "interaction zone re- 
gion" confined by this ring-like feature. If the in- 
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teraction zone region is filled with a denser, highly 

tTirbulcnt plasma, then it could potentially account 
for much of the dispersion towards PSR J1456-6843. 
If we assume an ambient density of » 0.01 cm~^ (i.e., 
a value comparable to that in quadrant 1), the in- 
teraction zone density {ug^iz ) will have to be ~40 
times larger in order to account for the excess dis- 
persion. Wc note that this is quite comparable to the 
value ~20-50(ne) derived from the SM estimate as 
constrained by our modeling. The consistency of the 
density estimates from the two independent methods 
can be argued, in some sense, to favor the existence 
of a dense interaction wall. 

Measurements of parallax for other nearby pul- 
sars in quadrant 4 will help to test for the existence 
of a wall. The most promising candidates for this 
seem to be low-DM (say, S 20 pc cm^'' ) objects with 
LOSs intercepting the interaction zone region. For in- 
stance, PSR J1751-4657 would be located much far- 
ther (-1765 pc) than its TC93 distance (-1080 pc), 
if the ambient Ug in this quadrant is comparable to 
that in quadrant 1 (~0.01 cm'^). PSR J1730-2304 
(close-by to PSR J1744 1134) is another interesting 
test case, whose distance will only be slightly more 
(«555 pc) if a denser zone is present, while it will 
be located much farther away (~1 kpc) if the entire 
sight line was uniformly filled at Hg ~0.01 cm~'^ . Sim- 
ilarly, the distance of PSR J1455-3330 will be ~950 
pc (compared to TC93 value of —740 pc) if we con- 
sider a dense zone and a low-density («0.01 cm~'^ ) 
ambient medium. While our arguments are based on 
a fairly simple picture, it is amply clear that measure- 
ments of a few interesting test cases will be valuable 
for understanding the LISM in this quadrant. 

3.3.4- Pulsars and the LISM: Further Prospects 

From the work presented here and related work 
in the recent past (see § 3.3.1), it is clear that pul- 
sars can be used quite effectively to probe large-scale 
features in the LISM. Conversely, it is also obvious 
that detailed interpretation of pulsar dispersion and 
scintillation data will need to take the structure of 
the LISM into account. There is growing observa- 
tional evidence that the Local Bubble is surrounded 
by numerous bubbles of similar properties, many of 
which are likely to be filled with hot X-ray-emitting 
gas. While this work has concentrated on Loop I, our 
most prominent such neighbor, other nearby exam- 
ples include the Eridanus bubble, the Gum Nebula 
and possibly radio Loops II and III. Further, there 
are also observations indicating possible interactions 
of some of these with the Local Bubble. In this paper, 
we have considered the possibility of an interaction re- 



gion between the Local and Loop I Bubbles. Toscano 

et al. (1999) and Chatterjce et al. (2000) have pro- 
vided some evidence for interaction of the Local Bub- 
ble with the GSH 238-^00-1-09 super-bubble in galac- 
tic quadrant 3. The closeness between the location 
of the near side of the Eridanus Bubble (159±16 pc; 
Quo et al. 1995) and the expected location of the Lo- 
cal Bubble shell in this direction (Ril30 pc) suggests a 
possible interaction between these two bubbles. The 
Vela supernova remnant has been shown to be embed- 
ded in a hot bubble confined by the shell of the Gum 
Nebula (Aschenbach, Egger & Triimper 1995). With 
the recently revised estimate for its distance (250±30 
pc by Cha, Sembach & Danks 1999), the near side of 
the "Gum-Bubble" is quite close to the Local Bubble 
boundary (wlOO pc); this is also supported by the ob- 
servation that the absorbing column density towards 
the Vela SNR is only about 10^" cm~^ . An interac- 
tion between these two bubbles therefore seems quite 
likely. 

Of the 163 pulsars known within ~2 kpc of the 
Sun, scintillation data are presently available for only 
73, and independent distance estimates for only 12. 
Clearly, a lot more needs to be done here. The new 
generation of large telescopes, such as the GMRT and 
the GBT, could potentially extend the available scin- 
tillation data for nearby weak pulsars. For many such 
cases, measuring pulse-broadening times at low fre- 
quencies (s 100 MHz) may prove to be a more vi- 
able technique than decorrelation bandwidth mea- 
surements. Observations of sources along well-chosen 
lines of sight where the scattering is dominated by the 
shells of these bubbles may also result in more detec- 
tions of unusual scattering effects, such as multiple 
imaging and extreme scattering events. Further, ISS 
studies of radio sources through cm- wave ISS and IDV 
may also provide useful insights into the nature of the 
LISM at higher Galactic latitudes. Many such obser- 
vations hold great promise for extending the rather 
simplistic models derived from the present investiga- 
tions into a more realistic picture of the LISM. 

4. Conclusions 

We have investigated the distribution of scatter- 
ing plasma in the LISM in the general direction of 
the Loop I Bubble by combining recent pulsar scintil- 
lation measurements. Many pulsars within ~2 kpc, 
and located towards 270° < I < S0° and -40° < /; < 
-1-80°, show enhanced levels of scattering, detected as 
a significantly reduced value of the ratio of the mea- 
sured to expected decorrelation bandwidths. The dis- 
crepancies cannot be explained by the Local Bubble 
model (BGR98) alone, and we interpret them as being 
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due to enhanced scattering associated with the shell of 
the Loop I Bubble and the Sagittarius spiral arm. Us- 
ing data from two heavily scattered pulsars with pre- 
cisely known distances, PSRs J1744-1134 and J1456- 
6843, wc have placed useful constraints on the scat- 
tering strength associated with the Loop 1 shell. The 
scattering measure inferred for the plasma inside this 
shell is ^0.3 pc m^^"/"^ , quite similar to that for the 
Local Bubble shell. Assuming a shell thickness ~1-10 
pc, this implies an average turbulence level ^100-500 
times larger than that in the ambient ISM. Adopt- 
ing this, our earlier model of the LISM is extended 
by incorporating an explicit scattering component for 
plasma in and around Loop I, in addition to that due 
to the Local Bubble. This two-bubble model suc- 
cessfully explains many of the observed decorrelation 
bandwidth discrepancies for nearby pulsars. The al- 
ternative possibility of enhanced scattering from an 
interaction zone between the Local Bubble and Loop 

1 is also considered. However, the level of turbu- 
lence in the interaction zone would have to be sev- 
eral times larger than that for the shells in the two- 
bubble model in order to explain the observations. 
Even then, the final agreement with observations for 
the interaction zone model is not as good as that for 
the two-bubble model. Assuming a simple relation 
between the scattering strength and the free electron 
density ( ocrig ) yields a shell density that is -^10- 
30 times larger than the ambient value, and a some- 
what higher density (~20-50(ne) ) for the interaction 
zone model. For several low-latitude pulsars at D^l- 

2 kpc {\z\:i 300 pc), we find that the observed scatter- 
ing discrepancies are consistent with additional en- 
hanced scattering from the Sagittarius spiral arm. 

We have discussed implications of our results for 
the interpretation of scintillation data as well as for 
the structure of the LISM. In the light of our results 
and those from several other recent works, we con- 
clude that the structure of the LISM needs to be con- 
sidered in the interpretation of pulsar dispersion and 
scintillation data, and may also be relevant for obser- 
vations of cm-wave ISS. Further, the general picture 
that emerges from the investigations in support of tur- 
bulent outer shells and/or interface regions of inter- 
stellar bubbles appears to be in qualitative agreement 
with the expectations of the models that describe the 
evolution of bubbles in the ISM. 
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A. Decorrelation bandwidth for an inhomogeneous, multi-component scattering medium 

Here we describe the method adopted for computing the decorrelation bandwidth {vd,pred)- Fig- 2 depicts 

typical scattering geometries relevant in our analysis: Db,LB is the expected location of the Local Bubble boundary, 
Dbi,Li and Db2,Li denote the nearer and farther boundaries, respectively, of the Loop I Bubble. D is the pulsar 
distance. 

For a homogeneous scattering medium, the decorrelation bandwidth [ud ) is given by (CWB85) 

2 

= ^ {Aa.f^.J^'^ C'^iz) dz^ (Al) 

where Aa is a model dependent constant, fobs denotes the frequency of observation, and a is the slope of the 
electron density wavenumber spectrum. In this scheme, the observer is at 2: = and the pulsar a,t z = D. If we 
adopt the canonical value of 11/3 for a (i.e., Kolmogorov-like spectrum), Ac = 2 x 10~^. 

The LISM model corresponding to Fig. 2 can be treated as an inhomogeneous scattering medium, with multiple 
components of different strengths of scattering ( , also called turbulent intensity or scattering strength), located 
at different points along the line of sight. The observable Vd is determined by the path length differences of scattered 
rays; therefore, contributions to it from C^{z) need to be appropriately weighted in such a way that scattering 
regions near the source or the observer produce smaller path-length differences than those that are mid-way 
(CWB85). Hence, for the case of our interest, equation (Al) can be rewritten as 



1 / 2 ^ I f^b,LB /•Db,LB +dLB 

l^d =J, {Aa^s)^^'' / M^) Cl^B^nMdz+ I W,{z) ClLBsh^z)dz+ 

r-Dbl.LI+C^L/ 




D 



Wci^z) Cl^i,^{z)dz+ I Wc{z) ClLish'y^)dz + 

Db,LB +dLB •^L'bl.LI 
/■Db2,LI /■Db2,LI+C?Z,/ 

/ Wo{z) Cl^LIint{z)dz+ Wc{z) Cl^LIsh{z)dz + 

"'l-'bl.LI+OL/ "^lJb2,LI 

I'D 1 (2^) 

/ , Cl^sm {Z)dz \ (A2) 

where dLB and d^i are the thickness of Local Bubble and Loop I shells, respectively. The turbulent intensity in the 
Local Bubble interior and in its shell arc represented by Cf^ lbuu ^^"^ LBsh ' respectively, and Limt 

L/s/t the equivalent quantities for the Loop I Bubble. The scattering plasma in the ambient ISM is charac- 
terized by . The symbol Wc{z) is the 'weighting function' for C^{z) , and is given by 

Mz) = i i^'i) (A3) 

This simple function is symmetric with respect to the mid-point between observer and pulsar, which means an 
inherent ambiguity is involved in the interpretation of the underlying scattering geometry. For instance, in the 
simple case of a thin screen placed between us and a pulsar, a screen closer to pulsar (say, at 2: = |Z)) will be 
equivalent to the one placed at z = jD from us. In the case of an extended, inhomogeneous medium (as in Fig. 2), 
the 'inverse geometry' is indistinguishable from the actual geometry. 
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Table 1: Pulsar Sample: The scintillation data 



PSR 


I 


b 


D 


^d^meas 


fobs 


Ref. 




[deg] 


[deg] 


[pc] 


[MHz] 


[MHz] 




(1) 


(2) 


(3j 


(4j 




(6j 


{') 


J1057-5226 


286.0 


6.6 


1530 


2.52 


660 


1 


J1430-6623 


312.7 


-5.4 


1800 


0.28 


660 


1 


J1455-3330 


330.7 


22.6 


740 


1.37 


436 


1 


J1456-6843* 


313.9 


-8.5 


455 


1.4 


436 


1 


J1537+1155 


19.8 


48.3 


680 


1.53 


436 


1 


J1543-0620 


0.6 


36.6 


1160 


0.111 


327 


2 


J1544-5308 


327.3 


1.3 


1290 


6.8 


1520 


1 


J1559-4438 


334.5 


6.4 


2000 


0.16 


660 


1 


J1603-7202 


316.6 


-14.5 


1640 


0.36 


660 


1 


J1605-5257 


329.7 


-0.5 


1240 


20.2 


1520 


1 


J1607-0032 


10.7 


35.5 


590 


0.379 


327 


2 


J1614+0737 


20.6 


38.2 


1500 


0.16 


436 


1 


J1709-1640 


5.8 


13.7 


1270 


0.040 


1000 


4 


J1713+0747 


28.8 


25.2 


1100 


1.45 


436 


1 


J1730-2304 


3.1 


6.0 


510 


0.17 


327 


3 


J1744-1134t 


14.8 


9.2 


357 


1.34 


436 


1 


J1751-4657 


345.0 


-10.2 


1080 


0.165 


327 


2 


J1752-2806 


1.5 


-1.0 


1530 


0.003 


1000 


4 


J1848-1952 


14.8 


-8.3 


960 


0.23 


436 


1 


J2053-7200 


321.9 


-35.0 


1110 


0.55 


436 


1 



References: (1) Johnston, Nicastro & Koribalski (1998), (2) Bhat, Rao & Gupta (1999b), (3) Gupta & Gothoskar 

(2000), (4) Cordes (1986) 

t Interferometric parallax distance from Bailes et al. (1990) 
^ Timing parallax distance from Toscano et al. (1999) 
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Table 2: Locations of the LB and LI boundaries and the IZ 



PSR 

(1) 






db.LB 

(2) 




dbi,Li* 
(3) 


db2,LI* 

(4) 




db,iz 
(5) 




J1057- 


-5226 


0, 


,04- 


-0, 


.05 


0.03 


0.14 





.03- 


-0 


.04 


J1430- 


■6623 


0, 


.03- 


■0, 


.04 


0.02 


0.15 





.02- 


-0 


.03 


J1455- 


-3330 


0, 


.05- 


-0, 


.10 


0.04 


0.42 





.04- 


-0, 


.07 


J1456- 


-6843 


0, 


.10- 


-0, 


.16 


0.07 


0.57 





.08- 


-0, 


.12 


J1537+1155 


0, 


.05- 


-0, 


.12 


0.07 


0.24 





.06- 


-0 


.10 


J1543- 


-0620 


0, 


.03- 


-0, 


.06 


0.03 


0.22 





.03- 


-0 


.05 


J1544- 


-5308 


0, 


.03- 


-0, 


.05 


0.02 


0.23 


0, 


.02- 


-0, 


.04 


J1559- 


-4438 


0, 


.02- 


-0, 


.03 


0.01 


0.15 





.01- 


-0 


.02 


J1603- 


-7202 


0, 


.03- 


0, 


.05 


0.02 


0.15 


0, 


.02- 


-0, 


.03 


J1605- 


5257 


0, 


.03- 


0, 


.05 


0.02 


0.24 


0, 


.03- 


-0, 


.04 


J1607- 


0032 


0, 


.05- 


0, 


.12 


0.06 


0.37 


0, 


.06- 


-0, 


.09 


J1614-F 


0737 


0, 


.02- 


0, 


.05 


0.03 


0.12 


0, 


.03- 


-0, 


.04 


J1709- 


1640 


0, 


.02- 


0, 


.05 


0.03 


0.19 


0, 


.02- 


-0, 


.04 


J17134 


-0747 


0, 


.02- 


0, 


.06 


0.06 


0.12 


0, 


.04- 


-0, 


.06 


J1730- 


-2304 


0, 


.05- 


0, 


.12 


0.07 


0.47 


0, 


.06- 


-0, 


.09 


J1744- 


-1134 


0, 


.07- 


0, 


.16 


0.12 


0.55 


0, 


.09- 


-0, 


.14 


J1751- 


4657 


0, 


.03- 


0, 


.06 


0.03 


0.24 


0, 


.03 


-0, 


.04 


J1752- 


2806 


0, 


.02- 


0, 


.04 


0.02 


0.16 


0, 


.02- 


-0, 


.03 


J1848- 


1952 


0, 


.03- 


-0, 


.06 


0.05 


0.16 


0, 


.04- 


-0, 


.06 


J2053- 


■7200 


0, 


.04- 


■0, 


.08 


0.05 


0.14 





.05- 


-0, 


.06 



Note: Columns (2) (5) list the fractional distances (e.g., db,LB =Db,LB /D). 

* For the Db,LB range of the solid and dashed envelopes in Fig. 4. 

* For the geometry in Fig. 4 (i.e., center located at ~170 pc from the Sun, towards (Z,6)=(329°, 17.5°)) 



Table 3: Various models considered and their goodness parameters 



Model Parameters eaii e<ifepc e>ifcpc 



LB 




SMlb =0.2 + ambient ISM 


1.06 


0.67 


1.41 


LB+LI 




SMlb =0.2, SMli =0.29 


0.55 


0.06 


0.90 


LB+LI-+-SA 


SMlb 


=0.2, SMli =0.29, C^^^ =0.047 


0.22 


0.06 


0.34 


LB-LI IZ-A 


SMiz 


=1.12, Db,iz =0.5(D^ LB +Dbi,Li) 


0.68 


0.13 


1.09 


LB-LI IZ-B 


SMiz 


=0.83, Dbjz =0.5(Db LB +Dbi,Li) 
SMiz =1.12, =0.047 


0.65 


0.11 


1.06 


IZ-A+SA 




0.27 


0.13 


0.38 


IZ-B+SA 




SMiz =0.83, C^'^^ =0.047 


0.25 


0.11 


0.37 
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Fig. 1. — Ratios of the measured decorrelation bandwidths (fd.meas ) to their predictions from the various models 

for the distribution of scattering in the Local ISM [vd.pred ) are plotted against the distance estimates. The 
Ud values are scaled to a common frequency of 327 MHz. The lone unfilled star indicates the measurement of PSR 
J1744-1134 at its TC93 distance of 166 pc. 
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Fig. 2. — Locations of the Local Bubble and Loop I shells along the sight lines towards PSRs J1744-134 and 

J1456 6843. The solid lines (Dbi.Li and Db2.Li) indicate the Loop I boundaries, and the dashed lines are the 
positions of the Local Bubble boundary (Db,LB ) that correspond to the two envelopes as shown in Fig. 4. The 
horizontal bar near the asterisk symbol (*) indicates the uncertainty in the distance estimate of the pulsar. 
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Fig. 3. — Same as in Fig. 1 for the model where the enhanced scattering is due to a possible interaction zone 
between the Local Bubble and the Loop I Bubble. The results are for Db,iz =0.5(Dg +Dbi,Li), i.e., for the 
IZ-A+SA model as described in Table 2 and S 3.2. 
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Fig. 4. - Geometry of the Local Bubble and the Loop I Bubble; panel (a) is the section in the Galactic plane, and 
panel (b) is the section along a plane perpendicular to the Galactic plane and passing through the north and the 
south Galactic poles, as well as through the centre of of Loop L The dashed and solid curves of the elongated cavity 
correspond to the two different geometries for the local scattering structure as derived from Ooty scintillation data 
(cf. BGR98). 
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